Ultrafast entangling gates between nuclear spins using photo-excited triplet states 
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The representation of information within the spins of electrons and nuclei has been powerful in the ongoing 
development of quantum computers 1 2 . Although nuclear spins are advantageous as quantum bits (qubits) due 
to their long coherence lifetimes (exceeding seconds 3 4 ), they exhibit very slow spin interactions and have weak 
polarisation. A coupled electron spin can be used to polarise the nuclear spirP* 7 -' and create fast single-qubit 
gates^, however, the permanent presence of electron spins is a source of nuclear decoherence. Here we show 
how a transient electron spin, arising from the optically excited triplet state of Ceo, can be used to hyperpolarise, 
manipulate and measure two nearby nuclear spins. Implementing a scheme which uses the spinor nature of the 
electron 10 , we performed an entangling gate in hundreds of nanoseconds: five orders of magnitude faster than 
the liquid-state J coupling. This approach can be widely applied to systems comprising an electron spin coupled 
to multiple nuclear spins, such as NV centres 11 , while the successful use of a transient electron spin motivates 
the design of new molecules able to exploit photo-excited triplet states. 



Different quantum systems possess different advantages as 
qubits, stimulating the use of so-called hybrid approaches 
to quantum computing 12 . Examples include interfacing su- 
perconducting qubits with spin ensembles 13 , optical pho- 
tons with defects in solidP, and electron spins with nuclear 
spins 4 ^^. Spins controlled with nuclear magnetic resonance 
(NMR) have played an important role in the developement of 
much experimental work in quantum information processing, 
showcasing high-fidelity control, complex demonstrations 
of quantum algorithms 17 and many-qubit decoupling strate- 
gies^. Nuclei in such systems are only weakly coupled: the 
indirect J-coupling interaction available in liquid-state NMR 
can be on the order of 100 Hz 17 , though nuclear spin dipole 
couplings in the solid state can exceed 10 kHz. 

This weak coupling places a lower limit on the duration of a 
quantum logic operation between two spins, and thus the com- 
putational speed of a nuclear spin-based quantum information 
processor. Additionally, the weak magnetic moment of nu- 
clear spins leads to a weak polarisation in general (typically 
less than 0.01% for liquid state NMR), making the scaling-up 
of the initial demonstrations very challenging unless a method 
for nuclear spin cooling can be applied 19 . 

These limitations can be addressed by making use of a cou- 
pled electron spin. Highly polarised electron spin states can be 
transferred to the nuclear spin coherently using SWAP opera- 
tion^Ellini or incoherently using a family of dynamic nuclear 
polarisation (DNP) methods^. Typical single-qubit gate 
times for electron spins are tens of nanoseconds, and given 
typical electron-nuclear couplings (in the range 1-100 MHz), 
it is possible to manipulate a nuclear spin on these timescales. 
For purely isotropic couplings, phase gates can be applied to 
nuclear spin qubits to perform dynamic decoupling 8 21 , while 
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using anisotropic coupling, more general gates have been ap- 
plied to single nuclear spins^"^^, or, recently, two nuclear 
spins^l. 

A disadvantage of using coupled electron spin is that the 
nuclear spin coherence time can be strongly limited by elec- 
tron spin relaxation or flip-flop processes 4 . A better strategy 
invokes an electron spin only at certain key times, for exam- 
ple to hyperpolarise the nuclear spins or to perform fast logic 
gates, so that there is minimal long-term impact on nuclear 
decoherence 2 ^. 

To explore such possibilities, we synthesised the fullerene 
derivative DMHFF 26 , illustrated in Figure [ij, containing two 
nuclear spins ( 31 P and 1 H) which are directly bonded to a 
Ceo fullerene cage. The molecule has a diamagnetic singlet 
ground state which can be photo-excited to populate the first 
excited single state. This state undergoes intersystem cross- 
ing (ISC) to a long-lived triplet state which is paramagnetic 
(S = 1) with electron spin density delocalised over the cage 
(see Figure [lj>,c). This system provides all the ingredients 
to explore nuclear spin manipulations mediated by a transient 
electron spin using a combination of optical excitation, elec- 
tron spin resonance (ESR) and NMR control. 

We begin by characterising the spin Hamiltonian of the 
DMHFP molecule: 

n=fi B Sg e B + SDS+ SA(i)Ii (1) 

*= 31 P, X H 

+ Jll.zh.z + li,nh ■ B 

where B the applied magnetic field, 7,4. „ the nuclear gyromag- 
netic ratio, g e the electron g e -factor tensor, /j,b the Bohr mag- 
neton, D the zero-field splitting (ZFS) tensor for the 5=1 
triplet state, A(i) the hyperfine coupling tensor between the 
triplet and the nuclear spins i, and J is the coupling between 
the nuclear spins. All terms involving S vanish in the elec- 
tronic ground state. 
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FIG. 1: The DMHFP molecule has an excited triplet state which 
couples to two nuclear spins: 1 H and 31 P. a) Illustration of the 
DMHFP molecule, b) Spin density distribution in the excited triplet 
state Ti, where blue areas correspond to atoms with high electron 
spin density, c) The system has a diamagnetic singlet ground state 
So which can be excited through ISC to a triplet state, Ti. The appli- 
cation of a magnetic field mixes the triplet sublevels in an orientation 
dependent manner. The triplet further couples to the two nuclear 
spins directly bonded to the cage, causing further splittings due to 
nuclear Zeeman and hyperfine energies, d) An electron spin echo- 
detected field sweep showing the ESR spectrum of the triplet state. 
Resonances at different fields correspond to different molecular ori- 
entations with respect to the applied field, e) Electron nuclear double 
resonance spectroscopy (ENDOR) is applied in four field positions in 
order to extract the isotropic hyperfine interaction between the triplet 
electron spin and the two nuclear spins 1 H and P, measured to be 
6 MHz and 1 1 MHz respectively). 



By performing pulsed ESR immediately following a 
532 nm laser pulse we examine the properties of the triplet 
state. Figure [TJi shows the intensity of an electron spin echo 
as a function of the applied magnetic field at X-band (9.7 GHz 
microwave frequency). By comparing the spectrum to sim- 
ulations^, we obtain the g e -tensor and the principal values 
of the ZFS tensor D xx = 52.13,£) TO = 159.53, D zz 

') 1 1 .fifi MHz. These ZFS parameters depend on the spatial 
distribution of the triplet wavefunction and characterise the 
strength and the asymmetry of the electron dipolar coupling. 

The hyperfine coupling between the triplet electron spin 
and the *H and 31 P nuclear spins can be measured us- 
ing the Davies electron-nuclear double resonance (ENDOR) 
method 20 , applied to select different molecular orientations 
within the sample. The spectra, shown in Figure show 
narrow peaks around 6 and 14 MHz corresponding to the nu- 



clear Larmor frequencies of the 31 P and 1 H spins, and aris- 
ing from nuclear transitions in the To subspace where the hy- 
perfine coupling is negligible (see Figure [lj). The ENDOR 
efficiency of these peaks is 90-100% enabling high-fidelity 
control on these transitions. 

The other peaks in the ENDOR spectra arise from the 
(orientation-dependent) hyperfine coupling in the T± sub- 
spaces. Fitting yields the isotropic hyperfine coupling terms 
A( X H) = 6.0 MHz and A( 31 P) = 1 1 MHz, consistent with DFT 
modelling (see Supporting Information). This hyperfine cou- 
pling allows for conditional electron/nuclear spin operations 
to be performed, however the breadth of these peaks (arising 
from the randomly oriented solid) results in poor fidelity nu- 
clear spin control in the T± subspaces. Nevertheless, we will 
show that is is possible to apply entangling operations in the 
To subspace where there is negligible coupling between the 
nuclear and electron spins. 

NMR experiments in the absence of optical excitation re- 
veal the ground-state J-coupling between 31 P and 1 H to be 
30 Hz, which leads to a nuclear controlled-NOT (CNOT) op- 
eration time of 17 ms. In the solid state the dipolar coupling 
can be measured using the spin echo double resonance (SE- 
DOR) pulse sequencers!], shown in Figure [2^. SEDOR can 
be interpreted as a standard NMR CNOT operation modified 
to allow for initialisation and readout by the electron spin. The 
timing of the refocussing pulses is swept to identify the CNOT 
coupling time of 160 /its, corresponding to a 3 kHz nuclear 
coupling. By comparing this coupling time to the lifetime of 
the triplet state (approximately 0.5 ms) and nuclear Ti times 
(0.20(4) and 1.9(4) ms for 1 H and 31 P respectively), we see 
that shorter entangling gate times are needed for higher fi- 
delity operations. This can be achieved by using a triplet elec- 
tron spin transition to apply an Aharanov-Anandan (AAj^®l 
controlled-phase (CPHASE) gate to the nuclear spins, as pro- 
posed in Ref 10. 

If a quantum state is taken through a closed-loop trajectory 
in Hilbert space, it acquires a geometric phase equal to half the 
solid angle mapped out by that trajectory. Because the field is 
on-resonance with an electronic transition corresponding to 
a particular configuration of the 1 H and 31 P nuclear spins, we 
can apply a Toffoli gate: a microwave pulse which only rotates 
the electron spin when the nuclear spins are, say, in the |ft) 
state. Applying a 27r pulse to the electron in this way imparts 
a 7T phase shift to the |fT) state, with respect to the others in 
the To subspace, equivalent to a CPHASE operation acting 
upon an effective nuclear spin basis. Both the CPHASE and 
CNOT operations are well-defined with respect to an effective 
eigenbasis where |4) is the nuclear spin eigenstate associated 
with the chosen electronic transition. The eigenstates |2) and 
1 3) differ from |4) by a *H and 31 P spin flip, respectively, and 
|1) differs from |4) by a flip of both nuclear spins. 

The duration of the CPHASE operation is limited only by 
the hyperfine coupling strength, such that a CPHASE gate on 
the timescale of hundreds of nanoseconds can be performed. 
To illustrate how this phase gate can be applied to the indi- 
vidual nuclear spins, we applied 27r microwave pulses while 
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FIG. 2: Two different ways to implement nuclear spin entan- 
gling gates using the triplet electron spin, a) A technique based on 
the spin echo double resonance (SEDOR) sequence can be used to 
measure the dipolar coupling between spins. The sequence resem- 
bles a Hahn-echo experiment on one spin ( H), however both spins 
are flipped during the refocussing pulse such that the sign of their 
coupling remains unchanged. Microwave pulses before and after are 
used to prepare and measure the 1 H nuclear spin coherence, b) The 
SEDOR sequence produces an oscillation in the nuclear coherence 
as r is swept, corresponding to a nuclear coupling of 3kHz. c) An al- 
ternative implementation of a CNOT gate consists of two Hadamard 
gates (7r/2 pulses) and a CPHASE gate created by applying a se- 
lective 27r pulse to the electron spin, d) An ultrafast nuclear phase 
gate (created by a 2n microwave pulse) is applied during nuclear 
spin Rabi oscillations on the 1 H spin at points marked with arrows, 
e) Comparison of the speed of the entangling operations. From the 
liquid state to the photo-excited solid state with the application of 
CPHASE gates there is an improvement of the entangling speed of 
five orders of magnitude. 



driving nuclear Rabi oscillations (see Figure |2jl). We veri- 
fied that this CPHASE behaviour was conditional by observ- 
ing uninterrupted Rabi oscillations on the complementary nu- 
clear subspace. A CNOT operation can be built by combining 
the CPHASE gate with Hadamard rotations (see Figure|2j;). 

In order to compare the performance of these two entan- 
gling CNOT operations, we attempt to put the two nuclear 
spins into a Bell state and then perform tomography of the ef- 
fective spin density matrix, building upon methods described 
elsewhere (see Supplementary Information). In short, each 
element of the density matrix must be mapped in turn onto the 
observable electron spin transition, through a combination of 
microwave and RF pulses. Notably, a CNOT (or similar) op- 
eration is needed when reading the zero- or double-quantum 
coherences, and these can be accomplished using either of the 
methods introduced above. To improve the fidelity of the to- 
mography, each element of the density matrix is imprinted 
with a particular time-varying phase applied to the nuclear 
spins using RF pulses. The total pulse sequence for gener- 
ating the pseudo-entangled state and measuring it using quan- 
tum state tomography is given in Figure [3^. 

In Figure [3J? and Figure [3J; we compare the density matri- 
ces obtained using the two different implementations of the 
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FIG. 3: Density matrix tomography results, a) General pulse se- 
quence for the extraction of the |1)(4|, |4)(lj elements of the den- 
sity matrix. The CNOT representations correspond to the entan- 
gling gates as represented by the shaded areas of Figure [2^ and Fig- 
ure [2|>. The pulses mapping a density matrix element to the observ- 
able changes depending upon the density matrix element, b) Density 
matrix obtained using the coupling-based entangling gates. The slow 
coupling leads to a fidelity of the operations of 34%. c) Density ma- 
trix obtained using CPHASE-based entangling gates. The ultrafast 
CPHASE gate entangles the nuclear spins in 220 ns and the fidelity 
increases to 65%. 



CNOT gate: 1) exploiting the nuclear dipolar coupling in the 
solid state, or 2) the triplet-mediated AA CPHASE operation 
combined with Hadamard gates. The fidelities of the final 
density matrices pu with respect to the ideal Bell state p%, 



calculated according to F(p B , p D ) = (Tr ^/pb P d ^/rb)) 
are 34% and 65% respectively. The increased fidelity of the 
latter approach is due primarily to the much shorter gate times: 
the CPHASE entangling gate is performed in only 220 ns, 
and adding the Hadamard gates yields a CNOT gate time of 
34.2 ps. In comparison, the CNOT based on the dipolar cou- 
pling had a duration of 160 ps. The maximum fidelity of each 
approach given the finite triplet recombination time for this 
molecule is 68% and 85%, respectively. The residual imper- 
fection is due to the limited fidelity of the CPHASE opera- 
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tion, as evidenced in Figure |2}i, and small gate imperfections, 
consistent with simulations incorporating a 4% error on each 
nuclear gate in the sequence. 

The lifetime of the double quantum coherence T 2 dqc ~ 
100 (is and the lifetime of the zero quantum coherence 
?2,zqc ~ 200 /is which indicates that they are limited by 
the hydrogen T 2 in the photo-excited state. Both nuclear co- 
herence lifetimes in the liquid state are longer by orders of 
magnitude, which motivates the further study of controlled de- 
excitation of the triplet state in place of the stochastic decay 
process to remove the electron as a source of nuclear decoher- 
ence. 

Although in disordered frozen solutions some limited ori- 
entation selection is possible, the use of a single crystal sam- 
ple would improve the uniformity and hence fidelity of the en- 
tangling operations. Due to the orientation-dependent triplet 
populations, a single-crystal sample would additionally allow 
for the selection of an orientation with the highest electron 
triplet polarisation, maximising the quality of entanglement. 
This technique is readily extensible to a wide range of molec- 
ular systems where a large number of nuclear spins are cou- 
pled to photoexcited triplet states, which can in turn be placed 
into arrays. 

METHODS SUMMARY 

Dimethyl [9-hydro(C 60 -I/ l )[5,6] fulleren-l(9H)-yl] phos- 
phonate was prepared following the procedure reported by 
Nakamura et al., Scheme 1^. Mono-functionalization of 
Ceo was performed using dimethyl phosphonate in a solution 
of toluene and HMPA at 120°C in the presence of oxygen. 
The product was purified by silica column chromatography 
(toluene and ramped to 10% ethyl acetate in toluene. 

Pulsed electron spin paramagnetic experiments were per- 
formed using an X-band (9-10 GHz) Bruker Elexsys680 spec- 
trometer equipped with a low-temperature helium-flow cryo- 
stat (Oxford CF935). The arbitrary phase radiofrequency 
pulses were generated using a using a Rohde and Schwarz 
AFQ100B together with an Amplifier Research 500W ampli- 
fier. Photo-excitation was achieved using an Nd-YAG laser 
at 532 nm with a 10 Hz repetition rate, using 10 mJ pulses, 
7ns in length. Microwave pulse lengths were 128 ns for tt/2 
pulses, and 220 ns for both it and 2tt pulses. The duration 
of RF pulses (both tt/2 and tt) was 17 /is. The samples were 
prepared in toluene-d8 with concentration of 4 x 10~ 4 M, de- 
oxygenated and flame sealed under vacuum and flash frozen 
in liquid nitrogen. 
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